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In  the  book,  "The  Special  Functions  and  Their  Approximations,"  a 
class  of  rational  approximations  for  the  generalized  hypergeometric  functions 
vas  developed.  Now  I^(z)  can  be  expressed  in  terms  of  a  qP-^  or  a  . 

Thus,  corresponding  to  each  form  and  a  choice  of  certain  free  parameters 
there  is  a  rational  approximation  for  I<j(z)  .  J.  C.  P.  Miller  has  shown 
that  Ini+v(^)»  ®  ^  positive  integer  or  zero,  can  be  approximated  by  use 
of  the  recursion  formula  for  I^+^Cz)  applied  in  the  backward  direction. 

If  this  scheme  is  used  together  with  each  of  two  certain  normalization 
relations,  then  rational  approximations  for  I^(z)  emerge  and  the",  rational 
approximations  are  identical  with  those  noted  above.  The  analysis  leads  to 
a  new  interpretation  of  the  baclr>fard  recursion  scheme.  We  also  study  a 
third  case  for  the  evaluation  of  Iittf<^(z)  ,  m  a  positive  integer,  by  the 
backward  recursion  process  which  presumes  that  I«^(z)  is  known.  In  each 
instance  a  closed  form  expression  for  the  truncation  error  is  developed 
which  leads  to  a  very  effective  a  priori  estimate  of  the  error.  For  each 
case  it  is  shewn  that  the  round-off  error  is  insignificant. 
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INTRODUCTION 


'  In  my  treatise  on  the  special  functions  ,  a  class  of  rational 
approximations  for  the  generalized  hypergecmetric  function  pFq  was 
developed.  These  approximations  depend  on  a  number  of  free  parameters. 

Since  I^(z)  caji  be  expressed  in  terms  of  a  or  a  ^^F^  ,  there  Is  a 

particular  rational  approximation  corresponding  to  each  of  these  hyper¬ 
geometric  forms  and  a  choice  of  the  aforementioned  free  parameters. 

^  '  <-■ 

The  idea  of  using  the  recursion  formula  for  I^(z)  in  the  back¬ 
ward  direction  to  generate  values  of  I^(z)  is  due  to  J.  C.  P.  Miller  [2]  . 

It  is  a  very  powerful  tool  and  the  notion  has  created  considerable  interest; 
see  [1,  Vol.  2,  pp.  159-166],  [3,4]  and  the  references  quoted  in  these 
sources.  The  Miller  scheme  together  with  two  certain  normalization  relations 
also  gives  rise  to  rational  approximations. 

In  a  conversation  Jerry  L.  Fields  conjectured  that  the  specific 
rational  approximations  noted  in  the  first  paragraph  ai'e  identical  to  the 
certain  rational  approximations  which  emerge  by  use  of  the  backward  recur¬ 
rence  scheme  noted,  i»4he  eeeceid  paragraph.  In  the  present  paper,  we,  verifj^^ 
this  conjecture.  In  addition,  we  develop  a  new  interpretation  of  the  Miller 
method,  Ue  also  study,  a  third  normalization  technique  which  is  sometimes 
used  with  the  backward  recursion  scheme.  A  ^osed  form  analytical  expression 
of  the  error  for  each  case  is  derived.  These  equations  are  valuable  as  they 
l.ead  to  simple  asymptotic  estimates  of  the  error  which  are  very  realistic 
and  easy  to  apply  in  practice.  It  is  demonstrated  that  the  round-off  errbr 
is  insignificant.  The  paper  closes  with  some  numerical  examples. 

In  the  main  body  of  the  ;i'aper,  we  find  it  convenient  to  deal  with 
the  modified  Bessel  function  Iy(z)  .  The  results  are  valid  for  all  z  in 
the  cut  complex  z-plane  -  n  <  arg  z  <  r  and  in  the  cut  complex  v-plane, 

|arg  V  I  <  n  .  In  this  connection,  we  should  note  that  I_^(z)  =  I^(z)  if 
V  is  an  integer  or  zero.  Thus  we  suppose  throughout  that  v  is  not  a 
negative  integer.  Actually,  it  is  sufficient  to  have  0  s  arg  z  &  tt/2  in 
view  of  the  definition  of  I.j(z)  .  Also  it  is  sufficient  to  have  R(v)  >  -1  , 
for  if  l2__.j(z)  and  -’’..^(z)  are  known,  computations  of  > 

m  =  a  positive  integer,  can  be  done  by  use  of  the  recursion  formula  for 
I-ia-vC^)  •  All  of  this  not  withstanding,  it  is  convenient  to  restate  some 
of  the  key  equations  to  facilitate  application  of  our  results  to  the  Bessel 
function  J.^(z)  .  This  is  done  near  the  end  of  the  paper. 
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RATIONAL  APPRCKIMATIONS  FOR  .Iy(z) 


Case  I. 

We  begin  with  the  representation 

(1) 

Tiieorem 

1. 

.  (z/2)''*n(2) 

Isjiz)  p(vfi)  h„(z)  ^  ®n(2)  - 

(2) 

■ 

n  (-n).  ('n+x)j^ 

♦n(z)  =  Z  (Z.U  ic>  •  > 

k=0 

(3) 

i 

1 

X  =  4/z2  ,  X  =  v+2-6  ,  6  =  0  or  6  =  1  , 

(4) 

3Fo(-n,n+X,l;X)  , 

(5) 

or 

lin(z)  =  (-)"n:(n+NH-2-6)nX"  o^(-2n-l+6-v;X'^)  . 

1 

(6) 

1 

H6re  for  convenience  we  introduce  the  notation 

i^f^'pizN  =  f- 

(7) 

and  our  usual  shorthand  notation  for  generalized  hypergecmetric  series 
applies,  see  [l,  Vol.^  1,  pp.  41,  42],  Further, 


I 


Sn(z) 


(z/2)^0(n^^) 

r(v+i)  ' 


=  1-6 -R(v)  if  -1  <  R(v)  s  1-6  ,  n  =  0  if  R(v)  i  1-6  ,  (8) 

I 

2 


/ 


■in(z)  =  (-)"nI(n+v<-2-6)nX"[l+0(n-l)]  ,  I 

whence 

v^)  “  •  '>('->  >  -1  ■  I 

and  so,  for  z  and  v  fixed,  R(\>)  >  -1,  the  approximation  process  is 
convergent . 


Proof:  Equations  (2)-(5)  follow  from  [l,  Vol.  2,  p.  96]  with 


a=0,f=g=0,p=0,q=l,p^=vfl, 
a  =  l-  6  )3-v,X  =  v+2-6  ,6  =  Oor6  =  l,Y~z 

O 

and  z  replaced  hy  z  /4  .  Notice  that  the  3FQ  series  for  hu(z)  in 
(5)  turned  around  is  the  alternative  form  for  hj^(z)  in  (6). 

Equation  (8)  follows  from  [l,  Vol.  2,  p.  lOS]  while  (9)  follows 
from  (6),  see  also  [l,  Vol.  1,  pp.  259-261],  and  (10)  is  now  obvious. 


Remark:  In  the  proof  developed  in  the  cited  source,  it  was  necessary  to 
suppose  that  R(v)  >  -1  .  Later,  we  present  a  new  formulation  of  the  error 
which  shows  that  v  is  unrestricted  save  that  v  is  not  a  negative  integer. 
So  throughout  this  work  v  is  arbitrary  except  as  just  indicated.  Compu¬ 
tation  wise,  the  exception  is  no  burden  since  I_n(z)  =  InC^)  • 


Theorem  2.  Both  tn(^)  satisfy  the  same  recurrence  formula 


hn(z)  +  (Ci+XDi)hn.i(z)  +  (C2+XD2)h„.2(2)  +  Csh^-sC^)  =  0  , 


,  -  (■gnn-2)(n-X-<-l) 

'1  (2n+X-4)(n+X-l)  ’ 


^2  (n+X-l)(n<-X-2)(2n+X-5)  ’  °2 


n(2ntX-2)(2n+X-l) 

n+X-1 


(n+X-l)(n+X-2) 


3 


(11) 


_  -(n-l)(n-.2)(2nn-2)(2nn-l) 

*^3  (2n+X-5)(2n+X-4)(n+X-l)(n+X-2)  ’ 


^rtiere  n  2  3  . 


Ptojj.  ;  See  [1,  Vol.  2,  Ch.  12]. 
Case  II.  Next  we  consider 


1n,(z) 


(z/2)V 

r(v+i) 


lFi(v+i;2v+l;-2z)  . 


Theorem  3. 


^  _  (z/2)^e2  =n(z)  _  V^) 

"  r(w+l)  gn(^)  ^  gn(")  ’ 

g  ^  (-n)k(n+2v+2)t^(v+^)^^  /.n+k,n+2v^-2+k,l 

=n(z)  2-  (2v+l)j^(v*-3/2)j^k:  3^1  \^v+3/2+k 

k"  0 


’22 


SnC^)  -  3^1 


/-n,n+2v+2,l 
l’'>(^v+3/2 


1' 
'2z]  ’ 


or 


n:(n+2Nr^2)  /-n-v-^  )  ) 

“  (2z)n(v<-3/2)  1^1  V-2n-2v-ir^/  ’ 


Further, 


_  (^2)  e!0i/l 

Wnvz;  - 


JU  -  l-2R(v)  if  -1  <  R(v)  ^  1  ,  (U  =  0  if  R(v)  2  h  • 


Also, 


en(z) 


nl(n^2vt2)ne^  (2^/4)  ^  (zg/32) (z^-8) 

(2z)''(v+3/2)^  r  ■  (n+w+l)^ 


0(n-^^) 


(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

,  (18) 
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whence 


Vn(2) 


(z/2)'’(2z)"(v+3/2)^0(n“) 

r(v+l)(n+2\»+2)j^n: 


R(v)  >  -1  , 


(19) 


and  for  z  and  v  fixed,  R(\))  >  -1  ,  the  approximation  process  is  con¬ 
vergent  . 


Proof:  Equations  (13)-(15)  follow  from  [l,  Vol.  2,  p.  96]  with 

a  =  0,  f  =  g  0  ,  p  =  q  =  1  ,  O'.,  =  =  2v+l  , 

Qf=l,3  =  2v,X  =  2v4-2  ,  y  =  z 

and  z  replaced  by  -2z  .  Equation  (16)  is  equation  (15)  turned  around. 
Equation  (17)  comes  from  [l,  Vol.  2,  p,  103]  while  (18)  comes  from  (16), 
see  also  [l,  Vol.  1,  pp.  133,  258-26l].  Thus  (19)  is  at  hand. 


Remark:  See  the  remark  after  Theorem  1. 

Theorem  4.  Both  3n(^)  Sn(^)  satisfy  the  same  recurrence  formula 

i 

gn(z)  +  (Ej^+2Fyz)g^_3_(z)  +  (E^+PF^/z  )g^_2(z )  +  E3g^_3(z)  =  0, 

_  -(n+v)  (n+2v-l)  _  4n(n+v) 

^1  (n+v-l)(n+2v'-l)  ’  ^1  n+2v+l  ’ 


-n(n-l)  _  -4(n-l)(n+v)(n+2v-l) 

^2  (n+2v) (n+2v+l)  ’  ^2  (n+2v) (n+2v+l) 

E  =  (n-l)(n-2)(n+v) 

3  (n+2v)(n+2v+l)(n+v-l)  ’ 


(20) 


where  n  a  3  . 


Proof:  See  [l,  Vol.  2,  Ch.  12]. 
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BACKWAED  RECURRENCE  SCHEMATA  FOR  GENERATING  Iv(z) 


The  technique  for  generating  I^(z)  by  use  of  the  recurrence 
formula  for  I^(z)  employed  in  tba  backward  direction  is  as  follows.  The 
recurrence  formula 


2(m+v+l)  ^  ^  ^ 


is  satisfied  by 


i(m+N))Tt 

Ini+v(2)  e  J^+s,(z)  > 


(21) 


(22) 


sin(m+v)TT 


-  Vv(^)J 


(23) 


In  this  work,  we  always  take  m  a  positive  integer  or  zero.  For  later  con¬ 
venience,  we  also  record  the  formula 


(z) 

L 

Let 

-  N 

q,(N)(z) 

^m,v 

^2  m*- v+3  I  ^m+2 , v  (m^■v+3)  'Pm+4,v'  '  *  '  ' 


(N)  ^ 

'Pn<-2,m  ^  ’  'P 


(N) 

N+l,v 


=  1  . 


(25) 


( 

Clearly  «?h  ®  linear  combination  of  the  solutions  (22)  subject  to 

the  conditions  (25)  and  we  readily  find  that 


+  e 


-in(N<-l-m) 


in  view  of  the  Wronskian  relation 


I^(z)K^l(z)  +  I^^(z)K^(z)  =  1/z  .  (27) 
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/ 


Suppose  that  we  are  given  the  normalization  relation 


03 

>(z)  =  H  • 


k=0 


Put 


(N)/  »  fN),  . 

lc=0 


and  consider 


C  (z)  = 

nr-v'  ' 


a(N)(z) 


,  m  s  IP-1 


We  can  now  prove 


Theorem  5.  lim  i^^^(z)  =  I  (z)  ,  0  s  m  <  lP-1 
-  m+v  m->-v 

N— >0D 


Proof:  Using  (26)  and  (29),  we  can  write 


,(N) 

im+M(z)  = 


9(z) 


N+m  l'f^v+2'  '  ni+v 


%+nP-2(^) 


ntv'  '  N+i  /z1 

E  E  (-)Wv(') 

k=0  *TP-v+2'^'  h=0 


(28) 


(29) 


(30) 


(31) 


(52) 


and  the  result  follows  frcan  the  known  behavior  of  the  Bessel  functions 
for  large  order.  That  is. 


mtv 


~  r(^'^l)  » 


(33) 


Vv(2)  =  i(z/2)““"''r(nt<-v)  [l+0(m"^)] 


(34) 


A 
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We  next  show  that  cp^*^^(z)  ,  m  =  N+1,N,N-1,.  .  .  can  be  repre¬ 
sented  in  terms  of  a  generalized  hypergecaietric  polynomial.  We  then  prove 
that  for  two  specific  choices  of  9(z)  ,  the  series  (29)  can  also  be  expressed 
in  terms  of  a  generalized  hypergeometric .polynomial;  and  further,  for  the 
two  choices  of  0(z),  respectively,  i^  ‘^(2)  and  the  rational  approxima¬ 
tions  and  ^j^(z)/gQ(z)  ,  respectively,  axe  equal.  Actually, 

we  first  state  and  prove  theorems  for  the  Case  I  situation  in  some  detail. 

The  corresponding  theorems  for  Case  II  are  stated  and  proofs  are  omitted 
as  the  details  are  much  aMn  to  their  Case  I  analogs. 

Another  choice  for  0(z)  previously  discussed  in  the  literature 
is  I^(z)  .  We  call  this  Case  III  even  though  the  corresponding  i^^(z)  = 

(z)  is  not  a  member  of  the  family  of  approximations  from 

which  Cases  I  and  II  were  derived.  We  defer  further  analysis  of  Case  III 
to  a  later  discussion  when  we  determine  closed  form  error  expressions  for 
all  the  cases. 


HYPERGEOMETRIC  REPRESENTATION  FOR  cp|^+.](z) 


Theorem  6. 


^N+ljV  N,v  z 

(z)  =  I .  ,  „(")  (z)  ■  Hr-'} . 

^N=l,v^  '  2  ^N-2,v'  '  z  ,3 

Z  2j 

(N)  ,  ^  12(rH-v-l)  16(lH-v-2) 

(N)  .  .  6(mv-l)  32(N+v-2)5  32(N<-v-3)5 

=  z  5  ’ 


(N) 


24(N+v-2)„  80(N+-v-5).  64(N4-v-4), 


'Pn-5,v^^)  =  2  ^ 

^  •» 


(N)  80(B+v-3)3  182(»*v-4)5  128(ttfv-5)7 

*  ,3  *  -- ;i —  *  —7 — 


etc. 


(35) 
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*  _ _ ^ 


ffi=0 


(2m+2n): 


2(k+l-6  )(2n-lrH+v)l^  / €6  -k,k+l4Tl-e6  ,k-2n-v-(l-6  )e  ,2n-k+2+v+e6  -6 


T  4 


^T\ 


V4, 


X  =  4/z^  ,  N  =  2n-6  ,6=Oor6=l,c=Oore=l, 


T1  =  e6  +  (l-e)(l-6)  ,  (36) 

where  6  =  0  or  1  according  as  N  is  even  or  odd,  respectively.  Also 

,  .  .  ,  .2n-ra+l-6  , 

J^hz)  =  (g/z)  r(2n+2-6+v) 

m,v  r(ni''v+l) 


m-1+6  irid6 
n+—^ —  ,-n+-^ 

.  -2n-l+6  -v,m+v+l,  -2n-l+6+m 


(37) 


where  [p]  is  the  largest  integer  s  p  . 


Proof:  By  induction:  The  Table  (35)  is  readily  developed  by  use  of  (21) 
and  the  starting  conditions  (25)  and  it  is  easily  verified  that  (36)  gives 
the  poljTicmiials  listed  in  (35).  Put  (36)  with  6  =  1  and  e  =  0  in  (24) 
(in  (24)  replace  m  by  2n-2k  )  .  Then  after  some  algebra,  it  is  seen 
that  the  coefficients  of  like  powers  of  z  vanish,  which  proves  (36)  for 
6=1  and  e  =  0  .  To  get  (36)  when  6  =  e  =  1  ,  use  (21).  The  case  6=0 
is  similar  and  we  omit  the  details.  Finally,  (37)  is  just  a  special  case 
of  (36).  To  connect  these  two  equations,  we  set  2n-2k+0  =  m  and  choose 
e  =  0  or  1  according  as  m  is  even  or  odd,  respectively. 


Remark  1:  We  have  given  more  polynomials  in  (35)  than  sui'e  necessary  for 
the  proof.  The  additional  entries  are  given  for  convenience. 


Remark  2:  If  in  (36),  e  =  6  =  1  ,  and  if  k 
and  v+1  ,  respectively,  then  we  get  (36)  with 
e  =  0  and  6  =  1,  and  if  v  is  replaced  by 
€=6  =  1. 
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and  V  Eire  replaced  by  k+1 
e  =  6  =  0  .  Again,  if  in  (36), 
v<-l  ,  then  we  get  (36)  with 


HYIERGECMETRIC  REPRESEMTATIONS  FOR  e^^^(z) 
AND  THE  FORMS  FOR 


Case  I.  Consider  the  normalization  relation 

oo  /  \k/ 


,  V  (-)^(2fc<-v)r(k+v) 

9(z)  =  ^(^;^-l)  ■  2—  r(v+i)k:  >  v  f  o  , 


k=0 


=  1  =  Iq(z)  4  2  21  (-)%k(z)  ,  N>  =  0  , 

k-1 

vhich  is  given  in  [l,  Vol.  2,  p.  45,  Eq.  (2)]  . 


(38) 


theorem  7. 


-  r(v.i)k:  »Jk7v<'> 

=  3F„(-n,nn,l;X) 


■  (2A)  - tt— > 


n: 


X  =  4/z^  ,  N  =  2n-6  ,  \  =  v42-5  ,  £  =  0  or  6  =  1  . 


(39) 


Proof;  We  consider  the  case  6=1  only  as  the  details  f t  r  6  =  0  are 
similar.  We  demonstrate  that  like  powers  of  X  in  the  suras  on  the  first 
two  lines  of  (39)  eu:e  equal.  Thus  we  must  show  that 


hj^(v) 


k 


Z 

r=0 


hr 


{~y  (2r4v)r(r+v)(-k)y(s+l+v)p 
rlr(v+l)(k!-l+v)j.(-s)j. 


(40) 


10 


/ 


> 


(-r(s-k);(vrt^i)^ 


n-m  =  k  ,  n+m  =  s 


The  case  k  =  0  is  trivial.  Assume  k  >  0  .  Clearly  Bjj(v)  is  a  poly¬ 
nomial  in  N)  of  degree  k  which  vanishes  if  v  =  -u  ,  u  =  l,2,..,,k  .  A 
straightforward  calculation  shows  that  =  0  whence  hjj(v)  also 

vanishes  if  u  =  l,2,...,k  .  Next  multiply  both  sides  of  (40)  by 
(k+v*-l)l^_l  .  Then  (40)  and  (41)  take  the  form 


h*(u) 


^  *  *  (-)^(2r^-v)^(r^-v)(-k)^(s+l+v)r^(2k^■v) 

■  Z_  ^r  »  ^r  r.T(v<-l)(-s)rr(k+l+v+r) 


X  (-)^(s-k):r(2x+v) 
~  s:r(\n'i) 


Now  each  side  of  (42)  is  a  polynomial  in  v  of  degree  2k-l  and  bJ^(v) 
vanishes  for  u=  l,2,...,2k-l  .  Also  hg(v)  vanishes  for  u=  l,2,...,k 
If  V  =  -k-v  ,  V  =  l,2,...,k-l  ,  then  bj  =  0  for  r  =  0,1,..., v-1  and 
b;^j  +  b^_j  =  0  for  j  =  0,1, . . .  ,(k-v)/2  .  So  h^(v)  and  Bj^(v)  have 
the  same  zeros.  Further,  it  is  easy  to  see  that  the  coefficients  of 
on  both  sides  of  (42)  are  equal.  Hence  h^(v)  = 
hjj(v)  =  B]^('j)  for  all  values  of  v  which  proves  the  theorem. 


Theorem  8. 


'pifv^z)  =  (2/z) 


1-6  {-)  (v+l)n+l-6 


♦n(2)  > 


N  =  2n-6  ,6  =  0  or  6  =  1  . 


Proof:  By  induction:  Using  (4)  aivd  (37),  we  can  readily  verify  the  state¬ 
ment  for  n  =  0,  1  and  2  .  A  straightforward  analysis  shows  that  both 
sides  of  (44)  satisiy  the  same  recurrence  formula  which  is  easily  deduced 
from  (11). 
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W 


From  Theorems  7  and  8,  and  (32)  and  (37),  we  have 


Theorem  9. 

(H)  tnt^) 

I  m-1+6  I  ^  m-l+fi  nr^-6 

(N)  (2/2 )  ^  ^ _ \-2n-l+6 -v),iit^\rt-l,-2n-l+6+m 

V.(-)  =  e^Cz)  r(.^.*LT 

(45) 

Thus  the  result  produced  by  use  of  the  recurrence  formula  for  l^{z) 
employed  in  the  backward  direction  together  with  the  nomalization  relation 
(38)  and  the  rational  approximation  given  in  (2)-(5)  are  identical 


Case  II.  From  [l,  Vol.  2,  p.  45,  Eq.  (5)],  we  have  the  normalization 
relation 


n(z) 


(z/2)''e^_  ^  (2k+2v)r(k^-2v)  ^  ^ 

r(v+i)  ^  r(2v+i)k:  ^k+v^^'’  t  ^  , 


=  =  Io(z)  +  2  X  Ik(z)  ,  V  =  0  . 

k=l 


(46) 


Here  to  avoid  confusion,  we  replace  0  by  n  in  the  notation  of  equations 
(28)  and  (29). 

Proofs  for  Theorems  10,  11  and  12  given  below  are  akin  to  those 
for  Theorems  7,  8,  and  9,  respectively,  and  we  skip  the  details. 


Theorem  10. 


N+1 


n(N)(z)  =  ^ 


k=0 


(2k+2v)r(k*-2v) 

r(2v+i)k: 


(z) 
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_  F,  /-N-l,If<-2NH-3,l 

“  (IH-1):  ^  ^  \s)+3/2 


(IH-I)I 


Su+i^^) 


I 


(47) 


Theorem  11. 


cpj^^^(z) 

0,N) 


(if^-1);  ®N+i^^^ 


(48) 


and  as  a  consequence  of  (44)  and  (48),  we  have 


(2\;+2) 


,1-6  ,  .n 


=  I!/!* - (-)_(v<-l)n.l.6  ,  (j)  ,  6  .  0  or  6  -  1 

(2n+l-6):  nJ  ’ 


Thecrem  12. 


(49) 


.(N),  ^  _  Qn^-i(^) 


m(N)/  > 

,  (H)r  .  ,  1  _  £ 

~  n^^^(z)  r(m+v+l) 


r  m-1+6 1  /  m-1+6  m+6 

,L  2  \(  ~^  2  ’"’^■^2 


ni+v  z  2  3 


,-2n-l+6+m,iTt<-v+l,-2n-l+6  -m 


F 

1  1 


2n+l-6  /  -2n-3/2-v+6 
V  -4n-3-2v+26 


2z 


(50) 


A  FURTHER  IHTERPRETATION  OF  THE 
BACKWARD  RECURRENCE  PROCESS 


write 


From  Watson's  treatise  on  Bessel  functions  [s,  p.  295],  we  can 


^2n+2-6+v^‘^)  ■■  ^2n-6  ,n>+2(2J)Iv^^)  "  ^2n+l-6  ,v+l^^)^v+l(2)  >  (51) 
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f  .  »  =  2"-®'  6  -  0  or  6  .  1  .  (52) 

where  'Pq^^(2)  is  defined  by  (37).  Consider  the  case  6  =  0  only.  We 
put  (51)  in  the  form 


Vl(^)  -  n(  )  -  I^^.i(z)T2n,^+2(2)  ’  ^  ■^2n,v+2(^) 

Further,  we  can  put 


(^) 


^2n+l,v-»-l^^^ 


E  (z) 
n 


^n(^)  "  F  (z)  * 


n' 


/  /  .Sn+'j+l 

_  ^2n+l,v+l 

r(2n+v*-2) 


(53) 


,  ,  ,2n+v+l 

_ ^2n,vt-2 

^n(2)  =  r(2n+-\)+2) 


Now  it  is  known  [s,  p.  302]  that  for  z  and  v  fixed,  z  /  0  , 


lim  En(z)  =  I^(z)  ,  lim  P„(z)  =  ^^(z)  .  (55) 

n  — *>flo 


From  (44)  smd  (52),  we  have 


Ejz)  -  Ln(2)  ,  Ln(z)  r.:(n+v+2) 


(56) 


idiere  )irn(z)  and  hj^(z)  are  given  by  (3)  and  (5),  respectively,  with 
6  =  0.  But  for  z  and  v  fixed,  z  f  Q  , 


lim  I^(z)  =  1  ,  (57) 

n— >00 


see  (9).  That  is. 


lim  Ejj(z) 
n— >00 


lim 
n— >oo 


(Ve)''t„(2) 

r(vti)V*) 


lira 

n->oi» 


(z/2)'q»J,”^(z) 

r(v^-l)9^”\z) 


Iv»(z)  . 


(58) 


11* 


/ 


A  similar  analysis  can  be  made  for  Fn(z)  .  Also,  a  like  study  can  be  done 
for  the  Case  II  scheme.  We  omit  the  details. 

Further,  the  above  shows  that  the  backward  recurrence  scheme  for 
the  computation  of  zly(z)/l^l(z)  is  the  same  as  the  well-known  truncated 
continued  fraction  representation  idxich  in  turn  is  the  same  as  the  main  , 
diagonal  Fade  approximation  for  this  function. 


ERROR  AHAiySES 


In  the  first  part  of  this  sectiop  we  develop  closed  form  repre¬ 
sentations  of  the  error  in  ij^^(z)  for  Cases  I-III  under  the  assumption 
of  exact  arithmetic.  This  type  of  error  arises  because  N  is  finite  and 
is  called  the  truncation  error.  From  each  analytical  representation  of 
the  error,  we  deduce  an  asymptotic  estimate  of  the  error  which  is  very 
realistic  and  easy  to  apply  in  practice.  The  results  for  Cases  I  and  II 
^eu  m  =  0  are  much  better  than  those  given  by  (10)  and  (19),  respectively. 
Further,  for  Cases  I  and  II,  if  z  and  v  are  fixed  and  n  is  sufficiently  ' 
large  with  respect  to  m  ,  the  relative  error  in  the  approximation  for 
Ijjtfy(z)  is  essentially  indei)endent  of  m-.  ' 

I 

An  analytical  formulation  of  the  round-off  error  is  developed  in  ' 
the  latter  part  of  this  section  where  it  is  shown  that  this  source  of  error 
is  insignificant.  i  i 

We  now  turn  to  a  study  of  the  truncation  errors. 


Case  I .  Let 


"  Vv(')  - 


(N) 


vrtiere  i,^^(z)  is  given  in  (45), 


(59) 

I 


( 

! 

'  '  ;  '  ' 

15 


r 


) 


Theorem  15.  If  v  is  not  a  positive  integer  or  zero, 

oF?(-2n-l+6-vjz^/4)Ei^^(z)  .  - r(yl-6+v)  j  p  /l  1 

0  1^  *  '  J  (n+l):r(2n+2-6+v)  '1  2yn+2,-tt<-6-v  I  '  / 


^ ^ ‘ ■  -— «-  -  -I  (z)l  (z) 

r(2n+2-6+v)  sin  vrr  -m-v'  ''  2n+2-6+v^  '* 

r(2n+2-6+v)  ^2n+2-6+v(^)^v(^) 


(-)  ttCz/2)  j 

■*■  r(2n+2-6+v)  sin  \»rT  ^2n+2-5+v^^^^m<-v^^^ 


(60) 


Equation  (60)  can  be  rearranged  so  that  with  the  aid  of  L'Hospital's  theorem, 
we  can  get  a  representation  of  the  error  when  v  beccanes  a  positive  integer 
or  zero.  We  do  not  give  this  result.  HoVrever,  for  arbitrary  v  ,  we  always 
have 


I 

/„\2n+2-6+v  /,  ,  .2n\ 

2(-)  (2/2) _  .  OUzA)^^^) _ 

r(2n+2-6+v)  l2n+2-6+v(2)^+v'^J  ■*■  r(2rd2-6+v)r(2n+-3-5+v)  ’ 


where  s  =  n-6+v(s=5os)  :  r  v  is  (is  not)  a  positive  integer  or  zero.  Clearly 

the  backward  recurrence  scheme  is  convergent.  Further,  for  n  sufficiently 
large,  n  »m  ,  the  relative  error  is  essentially  independent  of  m  .  For 
convenience  in  the  applicaiions  we  recordi  the  fonnula  ' 


(-)"(z/2)^”^^r('n^-l-6->-v)  ,  ,r  .  ,  .1,-1 

(rn-l):r(2n+2r6+v)  )J 


r(2n+2-6+w)  ^2nt2-6+w^  ' 


(62) 


i 
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Proof; 


We  have  need  for  the  formula  [l,  Vol.  1,  p.  2163 


la(z)lbfz) 


r(a+i)r(b+i) 


QFj^(a+ljz^/4)QF3^(b+l;z2/4) 


(z/2)^''''  /  (a+b+l)/2,(a+b+2)/2 

-  r(a+l)r(b+l)  2^3  a+l,b+l,a+b+l 


(63) 


where  it  mast  be  understood  that  none  of  the  numbers  a+1  ,  b^-1  ,  a+b+1 
is  a  negative  integer  or  zero.  Now  let 


r(m+v*'l)r(  -2n-l+6  -v) 
(2/2)-®"-^*^*” 


.  ,6  .  sj  ,  ,2n+2-6“m 

(-)  Trr(m-*-v+l)(z/2) _ 

r(2n+2-6+v)  sin  vrr 


^m+-v^^^^-2n-l+6-v^^^ 


w 

=  Z 

k=0 


2k 


Then 


(64) 


r  pv  T_1  /-k,2n+2-6+v-k  \ 

aj^  =  [2  ^:(-2n-l+6-v)k]  2^i(^m+v+l 

(  4-  Hli 

\"”  2  / 


m-H-6  ^ 

2  /k 


I  k 


( -2n-l+6  -v)j^(m+v+l)j^(  -2n-l+6  +m)j^ 


(65) 


whence 


ak  -  0 


m-1+6 


2n+l-6-m  , 


17 


®lc+-2n+2-6  -m  ~ 


24n+4-26-2iY^2n+2-6+v)r(2n+3-6+v) 


X 


/  m  .  3-6\  /  m  4-6\ 

-  2  2  j  -  2  ^  2 

(2n-t-3-6+N^)j^(  l-m.v)j^(2n+3-6 -m)  k; 


,  k  a  0 


(66) 


Thus 


v^)  =  24  ^  ^  ^ 


-2n-.l+6  -v,nt+-v^l,-2n-l+6+m 


^  \2nr^2— 6  •in. 


i-l^r.(m4-v+l)(2/2)‘ 

r(2n+2-6+v)  sin  vrr  ■^-m-v^2M2n+2-&+N)tzJ  >  (67) 


and  in  view  of  (23),  we  have  the  alternate  representations 


.  ,  ,  (-)^nr(m^-v+l)(z/2) 


2n+2-6 -m 


r(2n+2-6+v)  sin  \>rr 


Intfv,(2)  1^-  -  (-)  Kgn+l-S+v^^)  ^2n+l-6+v^^^J 


2^3 


-2n- 1+6  -  V  ,m+  v+ 1 ,  -2n- 1+6  +ni 
2n+2-6  -m 


^  (-)\r(ia^~’J+l)(z/2)*' 


r(2n+2-6+v)  sin  vrr 


^2n+2-6 +v^^ ^[tt  ’ 

(68) 


As  a  remark  aside,  the  combination  (26)  and  (68)  yields  (37)  and  so  we  have 
an  alternative  proof  of  (37).  Using  (45),  (59)  and  (67),  we  can  write 
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^F^(-2n-l+6-v;zV4)4”^2)  = 


.  (z/g) 


m+\> 


^(m+^/^■l) 


QF^(in+ vH  ;  z  ^4 )  qF^  ( -2n-l+6  -  V ;  z2/4  ) 


,  ,6  /  .,  /  ,2n+2-6-m 

(-)  _ T  (z)l 

r(2n+2-6+\))  sin  N>n  -m-v'  '  2n+2-6+v 


(z) 


]■ 


and  since 


oF5(-2n-l+6-v;z^/4)  -  QFi(-2n-l+6-v;z2/4) 


(-)’^r(n+l-6+Ni)  /l 
(n+l):r(2n+2-6+v)  1  2\n+2,-n+6 -v 


(69) 


we  readily  find  the  first  part  of  (60).  The  second  part  of  (60)  follows 
from  (23)  or  it  could  have  been  found  by  repeating  the  above  analysis  with 
the  second  equation  of  (68)  in  place  of  (67). 


Next  we  briefly  examine  the  situation  when  v  is  a  positive 
integer  or  zero.  With  v  =  r+e  >  the  ^.^2  "the  right-hemd  side  of  (60) 
can  be  expressed  as 


A" 

1  2\n<-2,- 


n+6  -V 


j^-6+r/l 

1  2  \n+2,-n+6 


-V 


S2M 


2k 


k=n-6+r+l 


(tt<-2)j^(-nt6  .v)j 


■6+ri 


r2 


\n+2,. 


•n+6  -V ' 


V.)- 


.  n+1+6 


(“)  TT(ntl)I(z/2) 


211+2-26 +2r 


sin  nvr(2n+3-6+r)r(n+l-6+v)r(l-e ) 


1 

^2n+3-6+r, 


1-e 


(70) 
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The  first  term  on  the  right-hand  side  of  (70)  is  defined  when  e  *  0  . 

When  the  second  term  on  the  right-hand  side  of  (70)  is  multiplied  by  the 
coefficient  of  the  2^F2  in  (60)  suid  the  result  is  combined  with  the  tern 
involving  Im+v(z)l2n+2-6+v(^^  (®0)>  i'*'  will  be  seen  that  one  can  pass 

to  the  limit  as  e— >0  .  The  final  expression  is  not  of  great  interest 
and  we  omit  further  details.  Equation  (61),  vdiich  is  important  for  practical 
considerations  (see  later  numerical  example)  readily  follows  from  (60)  and 
the  above  remarks,  and  (62)  is  but  a  simplified  version  of  (61), 


Remark;  Let  v  ,  n  and  z  be  fixed  so  that  ^ui^ction  of 
m  only.  Then  E^^^(z)  satisfies  the  recurrence  fomula  for  ^(z),  see 
(21).  This  is  evident  from  (30)  and  confirmed  by  (60). 


Case  II .  Let 


,(») 

m,v 


(z) 


I 

m+v 


(z) 


(K). 


i'7(z) 

m+v 


irtiere  i^y(z)  is  given  in  (50). 


(71) 


Theorem  14.  If  neither  v  nor  is  a  positive  integer  or  zero,  then 


-z  2n+l-6 [ -2n-3/2+6 -v 
®  1^1  \.4n-3+26-2v 


+ 


(2z)2n+2-6(%^i) 


-z 


n+2-6 


(2n+2-6):(2n+2-6+2v) 


2n^■2-6 


2^2' 


i“V,l 

2n+3-6 ,-2n-l+6-2v 


r(2n+2-6+v)sin  vrr 


(72) 


idiich  is  the  same  as  the  right-hand  side  of  (60). 

If  is  a  positive  integer  or  zero,  call  it  r  ,  we  have 
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:  2n^l-6/-2n-2+5-r 

ri  \-4n-4+26-2rl  /  m+r+a' 


(r+l)2n+2-6 


(2n+2-6  ):(2n+3-6+2r) 


2n+2-6 


®  Vr+|^^^2^(2n+3^ 


.2n-2+6  -2r 


,  fi+r+1  ,  ,  2n+5/2-6+r^ 

(-)  tt(z/2)  I, 

r(2n+5/2-6+r) 


2n+r<-5/2-6 


r  =  0,1,2,...  ,  4  =  1  if  r  =  -1  ,  p,  =  2  if  r  2  0 


In  particular, 


(-}^  (tt/2z)%~\ 


2n<-3/2-6' 


^n<-3/2-6^^^ 


,  .6,  .2  -z,  _  ,  _ 

=  (-)  (grrz)  e  (z/2) -  [l+0(n‘^)]  .  (74) 

[(2n+l-6):] 

As  in  Theorem  13,  (72)  can  be  rearranged  to  get  a  representation  for  the 
error  when  v  is  a  positive  integer  or  zero.  This  result  is  knitted • 
However,  for  arbitraiy  v  ,  v  /  -i  >  we  always  have 


2ntl-6(-2n-3/2«.>,L  \  (H) 

'  1^1  \.4n-3+26-2Ml^^/*in,v‘^^ 


,  2n+2-6  .  1 , 

2n^-2-6  ^-z-r  (z)  Uz 

(2n<-2-6):(2n+2-a+2v)2^^2_6  2  2\2n+3-6  ,-2n-l+6-2vl  j 


,  m*-6  .  ,  2n+2-6+v 

*  ■  ’  - W-s*v(*)VvU) 


,  ou^/g)  )  , 

r(2n+2-6+v)r(2iH'?-6+v) 
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^ere  t  =  \>-2(t=2n+l-6+2v)  if  v  is  half  a  positive  integer  (is  a  posi 
tive  integer  or  zero)  and  where  t  =  »  for  all  other  v  .  Clearly,  the 
cackward  recurrence  scheme  is  convergent. 


Further,  for  n  sufficiently  large,  n  »  m  ,  the  relative  error  is 
independent  of  m  .  For  convenience  in  the  applications,  we  record  the 
formula 


f(n) 

m,'j 


(z) 


_  ,2n+2-6,  1.  -z 

-(2z)  ('^■^2)2n+2_6® 


(2n+2-6  ) :  (2n+2-6+2v)2jj^2_6 


r(2n+2-6+v) 


Proof :  Let 


(-2n-3/2^6-v|  \1.  ,  f- 

ri\-4n-3>25-e«'  /  f— 

J  k=0 


hj^^Z 


(77) 


In  view  of  (1)  and  (12),  Bj^(z)  =  \(z)  \rtiere  Anj(z)  is  given  by  (64). 
Hence  b^  =  0  for  k  odd  and  b2k  =  ajj  where  aj^  is  given  by  (65)-(66). 
The  analysis  proceeds  as  for  Case  I  and  we  find  (72).  Notice  that  the 
right-hand  sides  of  (72)  and  (60)  are  identical. 

For  the  proof  of  (73),  let  Vj„  ^(z)  be  the  entire  first  term  on 
the  right-hand  side  of  (72),  that  is,  the  tern  involving  the  2^2  •  Le'^ 

V  =  r+-^e  ,  r  =  0,1,2,...  .  Then  we  can  write 


=-Cv(z)(r+l+€)2n+2-e 


f  rrfi-''-! 

0*0'  On+'2 


2nt2-6+2r 

*  Z 

k=r+l 


2'2\2n+3-6  ,-2n-l+6  -2v 
2k  1 


2z 


2k  ^ 

CkZ  +  2_ 

k=2n+3-6+2r 


Cj^Z 


C.,(z)  = 


(2z)' 


‘e-^I  .  (z) 
m-’V 


(i-v) 


(2n+2-6):(2n+2-6+2v)2jj+2_6  ’  (2n+3-6  )j^( -2n-H-6 -2v) 


It  is  easy  to  pass  to  the  limit  when  e — >0  and  (73)  residily  follows. 
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Now  for  arbitrary  a  and  b  , 


(2/rrz)^(a  cosh  z  +  b  sinh  z)  , 


Further,  from  (16),  with  n  replaced  by  2n+l-6  ,  see  also  (17),  ajid 
V  =  -^  ,  gn(z)  reduces  to 


^o(- 


2n-l+6,2rr^2-6;-— \=  (rr/2z) 
2z  / 


®  ^2n+3/2-6^^^ 


Use  these  data  with  r  =  -1  and  m  =  0  in  (73)  to  get  the  first  line  of 
(74).  Derivation  of  the  second  line  of  (74)  is  trivial  and  details  are 
omitted.  Since 


I_l(z)  =  (2rTz)~^e^(l+e"^^)  , 

-  -z  (N) 

it  is  easy  to  show  that  (2TTz)^e  i  i,  (z)  is  the  main  diagonal  Fade  approxima- 

tion  to  1  +  e"  •  When  allowance  is  made  for  a  change  of  notation,  (74) 
is  a  previously  obtained  result  [l,  Vol.  2,  p.  74,  Eqs.  (34),  (35)]. 

When  V  =  i*+-e  ,  r  a  positive  integer  or  zero,  we  can  rearrange 
(72)  after  the  manner  of  the  discussion  surrounding  (70)  and  use  L'Hospital's 
theorem  to  get  the  limit  as  e  — >0  .  The  result  is  not  of  immediate  interest 
emd  we  omit  detad-ls.  The  statements  (75)  and  (76)  are  readily  derived  and 
here  too  we  skip  details. 


(  n) 

Remeurk:  F'  '(z)  with  n  ,  v  and  z  fixed  and  m  variable  satisfies  the 
—  m,v 

recurrence  formula  for  ^(z)  ,  see  (21).  This  feature  is  clearly  depicted 
by  (72).  In  both  (62)  and  (76),  the  term  involving  K^^(z)  is  of  lower 
order  than  the  term  involving  ljjj+y(z)  •  Neglecting  the  former  term  in  each 
equation,  we  have 
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Ilheorem  15. 


m,v 
m,v'  ’ 


(.)“"V(n+l-6+v)(2nf2-6  ):2®^‘^(i)^(z/2)^* 
(n+l)lr(2n+2-6+2v) 


.  vH+l, 1 .  .  .  .8  z 

(-)  (2L(z/2)  e 

-  [l+0(n-^)]  ,  v  ff  -i  . 


n 


v*-6 


(78) 


This  shows  that  there  is  little  difference  in  the  accuracy  of  the 
two  schemes  for  the  evaluation  of  Intt.,j(z)  .  Computation-wise,  if  the 

hackward  recursion  scheme  is  used,  Case  I  requires  less  operations  since 
the  associated  normalization  relation,  see  (38)  and  (39),  uses  the  sequence 

)}  ,  k  =  0,2,4,...  ,  ^rtiile  the  Case  II  normalization  relation,  see 

(46)  and  (47),  employs  ^  ~  0,1,2,...  .  Also  to  get  I^(z) 

by  the  Ceise  II  scheme,  e^  must  be  evaluated.  On  the  other  hand,  if  |z| 
is  large,  r(z)  >  0  ,  one  often  wants  not  I^(z)  ,  but  e”^I^(z)  .  The 
latter  is  automatically  funiished  by  the  Case  II  technique.  It  appears  that 
for  the  same  n  ,  the  Case  II  procedure  might  be  more  accurate  than  the 
Case  I  scheme  even  for  moderate  values  of  |z|  ,  R(z)  >  0  ,  in  view  of 
the  presence  of  e^  in  the  numerator  of  (78).  Also,  Case  II  is  favored 
when  R(v+6 )  <  0  .  Improved  information  cannot  be  derived  from  (78)  as 
the  estimate  is  for  fixed  m  ,  v  and  z  .  For  error  analyses  it  is  sug¬ 
gested  that  one  use  (62)  or  (76)  as  appropriate.  Farther  discussion  is 
deferred  to  a  later  part  of  the  paper  where  numerical  examples  are  presented. 

If  z  is  pure  imaginary  and  v  is  real,  then  z  '^I^(z)  is  real 
and  definitely  the  Case  I  procedure  is  better  than  the  Case  II  scheme  since 
the  former  requires  real  arithmetic  while  the  latter  demands  complex  arith¬ 
metic  . 


If  only  Iy(z)  or  only 


e-"l. 


,(z)  is  required,  use  of  the 


rational  approximation  scheme  or  the  equivalent  backward  recursion  scheme 
demands  about  the  same  number  of  operations.  In  the  absence  of  a  priori 
estimates  of  the  error,  the  rational  approximation  scheme  employed  in  the 
following  fashion  is  preferred.  It  is  sufficjont  to  consider  the  Case  I 
situation.  Compute  ^^(z)  from  either  (3)  or  the  combination  (37),  (44), 
and  Mz)  from  (5),  for  n  =  0,1  and  2  .  Compute  subsequent  values 
of  tn(z)  a»i  ^^(z)  by  use  of  the  recursion  formula  (12).  Comparison  of 
f^(z)/hjj(z)  with  ^^^(z)/h^_^^(z)  affords  an  estimate  of  the  error.  If 


rin-iv“//“n+J' 

one  requires  Ij^^.^(z)  or  e  Ij^y(z)  for  k=  0,1,2,., 
the  backwai’d  recursion  scheme  is  highly  advantageous. 


,r  ,  then  obviously 


Case  III .  Consider 
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Finally, 


for  convenience  in  the  applications,  we  record  the  formulas 


m,v 


(z) 


m+S  2n+4-26+2v-m 
(-)  (z/2)  r(m+v) 

r(2n+2-6+v)r(2n+3-6+v)r(v+l) 


F 

2^3 


Wm-1  m-2 
2  ’  2 
-v-m,\>+l,l-m 


[l+0(n-l)] 


1^(2) 


(83) 

(84) 


Proof:  Using  (37),  (67)  and  (69),  we  find 


cp^  ^^(z)G^  ^(z)  =  ''  /.  ^2n+2-5+v(^^  ' 

UjV  DijV  Sin  vrr 


smd  in  view  of  (23) 


’  <-)*  W2-6*s,<^>  {'  • 


Frcm  (26),  the  portion  in  curly  brackets  in  the  latter  equation  is 
tn  {m  2  ^ 

(-)  ©'  "  ■'(z)  ,  Hence  the  first  line  of  (81)  is  at  hand.  The  remainder 

'  ^0,v 


of  (81)  follows  from  (37).  The  first  line  of  (81)  coupled  with  the  dis¬ 
cussion  surrounding  equations  (21)-(25)  produces  (82).  By  the  confluence 
principle,  see  [l,  Vol.  1,  p.  50] 


,  ,  N+1  ,  ,  fN+1]  /  N+1  N 

(N),,,  ^  ^  Iti  A— 


cp'”'(z) 


r(v+i) 


2^3  \ 

\_N-1-v,v+1,-N-1 


N^l 


=  (2/2)^^'^\(lf^v+2)l^(z)[l+0(r^)J  , 


(85) 


/ 


t 


whence  (83)  and  (84)  readily  follow  with  the  aid  of  (33). 


Error-wise,  it  is  difficult  to  compare  Cases  I  or  II  with  Case  III 
without  some  simplifying  assumptions.  If  n»  m  ,  using  (33),  (62)  and 
(83),  we  get 


,2m+2-26+2\j-m 


^  (-)  (z/g) _ 

„(N).  .  r(2n+3-6+v)r(n+l-6+v)r(v+l) 


(n~*-l)ir(m+\;) 


fm-1 

/  m-1  m-2 

L  2 

f  ■  2  2 

’’3 

\l-v-m,sH-l,l-m 

.[l+0(n-^)] 


(z/2) 


2n 


r(2n+3-6+v) 


6+1-v 
0(n  ) 


(86) 


and  so  Case  III  is  superior  to  Case  I.  Now  suppose  m  is  sufficiently 
large  so  that  in  (62),  the  second  term  dominates  the  first  term.  This  is 
certainly  the  case  if  m  =  2n+i-6-d  ,  d  «  n  ,  in  view  of  (33)  and  (34). 
Then 


m,v^  ' 
m,v'  ' 


=  (2/z)®  [l+0(n"^)l 


(87) 


and  under  these  conditions  there  is  little  to  choose  between  the  two  cases. 
Overall,  it  appears  that  Case  III  gives  better  accuracy  than  Case  I.  How¬ 
ever,  for  Case  III,  one  must  know  Iiy,(z)  vdiile  for  Case  I  no  such  knowledge 
is  required.  For  all  z  >  C  and  all  v  ,  0  s  v  s  1  ,  coefficients  are 
available  to  facilitate  the  rapid  evaluation  of  Jy(z)  and  I^(z)  ,  see 
Ql,  Vol.  2;  6,7,8j.  (Actually,  much  more  information  is  given  in  these 
sources.)  All  of  this  csin  often  make  the  Case  III  approach  rather  attrac¬ 
tive.  See  the  numerical  examples. 

Next  we  consider  the  round-off  error. 


Cases  I, II.  It  is  sufficient  to  trace  the  effect  of  a  given  round-off  error 
in  a  single  entry  of  the  table  generated  by  the  backward  recursion  process. 
Thus  let  ZVp  be  the  symbol  for  the  round-off  error  in  cp  .  Suppose  that 
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A«P^^^(z)  =  0  for  m  =  N+2,N+1,...,  S+2  , 
m,v 


=  (ju  for  m  =  S+1  ,  (88) 

where  S  =  2s-y  and  Y  is  0  (is  1)  if  S  is  even  (is  odd).  Hence 


(89) 


Ae^^\z)  =  (uG^^^z)  . 


(90) 


For  the  Case  I  procedure,  we  have 

0(N)(^)  ■  e^N)(z)  -U)e^®\z) 

'j  r  “’oW, 


{e"''(^)l 


e^^\z) 
0W(z) 


(juoE^^\z) 
^  m,v'  ' 


E^C.) 


jl  -  m,v 


E^^\z) 


IttjV 


l-UjO 


(91) 


and  these  equations  also  hold  for  Case  II  provided  E  is  replaced  by  F  . 
If  S  =  N  ,  that  is,  s  =  n  ,  the. round-off  error  is  nil.  Indeed,  this  must 
be  since  the  starting  value  cp  is  immaterial.  It  is  clear  that  if 

an  parameters  and  z  are  fixed,  then  the  round-off  error  decays  to  zero 
as  n— >»  .  For  n  sufficiently  large  with  respect  to  z  ,  (a|<  1  and 
we  can  take  |l-E^^^(z)/E^®^(z)j <  2  .  Also  it  appears  heuristically  that 


E^^\z)  <|e^^^(z)  s  It  (z)  -  (z/2)7r(v+l) 

m,v  0,v  \j 


(92) 


28 


Thus  on  this  basis, 


AVy(z)l^  -  (z/2)Vr(v+i)| 


(93) 


is  an  approximate  hound  for  the  round-off  in  a  single  entry  of  the  set  of 

numbers  generated  by  the  backward  recursion  process.  If  tu  is  the  m^imum 

round-off  error  in  each  entry,  then  the  total  round-off  error  in  i'^Mz)  is 

m+v 

approximately  bounded  by  N  times  the  right-hand  side  of  (93),  Thus  the 
round-off  error  is  insignificant,  and  it  is  easy  to  estimate  the  numt^er  of 
extra  decimals  which  must  be  carried  so  that  the  total  round-off  error  in 
the  process  lies  within  the  error  when  the  arithmetic  is  exact.  Equations 
analogous  to  (92)  and  (93)  for  Case  II  are  easily  derived  and  we  omit  details 


Case  III.  We  have 


^  m,v^ 


ra,v  0,v  m,v  0,v  m,v 


Using  (26),  a  straightforward  computation  shows  that 


(94) 


N,s  y  2r  6+y+i  ■) 

[^2s+2-y+v  ( 2 )  ^n+2-6  +v(  ^  )  ^2n+2-6  +v  ( ^  )^s+2-Y+v  (^  )J 


m+Y+1  (N)  (m-2) 

‘(’2s+2-y,v'P0,v  (z).'">0 


=  0  ,  m  =  0  . 


(95) 
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’o.v  ''ni.s.'  (»)(,) 

D,v 


Obviously,  this  is  nil  is  S  =  N  or  if  m  =  0  .  .  Now  using  (37),  we  find 


\"/i  f_\  ,  ,  2n+2s+3-6 -Y-m  ,  .  ,  . 

^2s+2-v,v^^^  ^  (2/2) _ r(v+l)r(m+v+l) 

cp(”^('z)cp(^^(z)  ^  r(2n+2-6+M)r(2s+3-6+v) 


0,v  ’  m,v 


n-s- 


F 

2  3  - 


2  J  2  ’  2 

\-2n-l+6  -v  ,2s+3-Y+NJ ,  *-2n+2s+l+6  -Y 


6-11/  6-1  6 


X  ra-1+61  /  m-1+6  nFs 

2\  v[«-“T“/-^^“2"’-"^ 


-n-— ,-n-^  z^\  fL  2  j/  ‘‘  2  ’  “  2 

-2n-l+6-v,v+l,,-2n-l+6  \-2n-l+6 -v,m+v+l,-2n-l+6+m| 


and  it  is  clear  that  rounding  errors  are  insignificant.  Indeed,  if  all 
parameters  and  z  are  fixed,  then  the  round-off  error  tends  to  zero  as 
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FORMULAS  FOR  J^(z) 


As  previously  remarked,  the  analyses  for  I^(z)  hold  throughout 
the  cut  ccanplex  z-plane,  -n  <  arg  z  s  rt  ,  and  throughout  the  cdt  ccmplex 
v-plane  (arg  v(<  tt  >  although  it  Is  sulTicient  to  have  0  s  arg  z  s  tt/2  and 
R(v)  >  '1  .  Nonetheless,  ve  indicate  how  to  get  results  for  Jv(z)  directly 
eind  to  facilitate  use  of  our  findings,  it  is  convenient  to  restate  some  of  the 
key  equations.  We  omit  discussion  of  Case  II  since  it  requires  ccanplex  arith¬ 
metic  to  generate  Jy(z)  which  is  real  when  z  and  ‘j  are  real.  In  einy 
event,  the  reader  should  have  no  difficulty  in  establishing  the  Case  II  equa¬ 
tions  for  J.^,(z)  once  it  is  observed  how  this  is  done  for  Case  I. 

All  developments  for  J^(z)  are  readily  gotten  by  use  of  the 

equations 


-i(i3<-v)n/2 


.  ,  -W2,  ,  (1) 


m+v' 


Y^^(z)  =  (-)“csc  vrr[(-)°cos  vn  Jnifv(2)‘J-m-v(2)]  » 


(96) 

(99) 

(100) 

(101) 


where  now  in  the  J^(z)  analyses,  -tt/2  <  arg  z  s  5n/2  . 

It  is  convenient  to  introduce  the  following  notation.  CMless 
indicated  otherwise,  if  A  is  used  to  signify  some  function  or  equation 
in  the  developments  for  I.^(z)  ,  then  A*  is  used  to  signify  the  correspond¬ 
ing  function  or  equation  in  the  developments  for  Jv(z)  .  In  illustration 

J^(z)  =  (z/2)  QFifv+l;-z^/4)  ,  (1)* 

'and  both  Jnjfv(z)  and  Yjj+^(z)  are  solutions  of  the  difference  equation 


*  ,  .  2(m+v+l)  *  ,  X  * 


(21)* 
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Also 


,in+v 


Jr:L!A>(2)  =  ffavl+i]  l+0(m-^)  , 


(33)^ 


'  -  „ —  r(m»v)  [1+0(11.-^)] 


m+v 


(34)* 


We  now  present  the  key  results  pertinent  to  J^(z)  . 


Theorem  1*. 


(z/2)^ 

Jsj(z)  =  p(m+i)  *  /  Rn(^)  »  Rn(^)  ~  ~ 


S„(z) 


h  (z) 
n' 


(2)* 


where  tn(^)  hj;(z)  are  given  by  i^q(z)  and  h^Cz)  ,  see  (3)-(5), 

provided  there  we  replace  z^  by  -z^  ,  that  is,  replace  X  by  -X  . 
Further,  S^(z)  and  8^(2)  are  both  given  by  (8)  and  Rn(z)  =  (-)”Rn(z)  , 
see  (10). 


Theorem  2.*  Both  i|tn(z)  and  hji(z)  satisfy  the  same  recurrence  formula 
(3J.)  if  there  we  replace  X  by  -X  . 


We  state  without  proof  the  following  equations. 


(W)* 


(z) 


( z/2  p  ( 2n*2-6  +v  ) 

r(ni+v+l) 


-n+ 


u*— A  •  y 


,-n+ 


2  2 
S^-2n-l-  v*5  ,in+v+l,-2n-l'^6' 


(37)* 


e*(z) 


r(v+i) 


=  f- 

^  r(v+i)k: 
k=0  '  ’ 


(38)* 


e 


(2/z) 


1-6 


iliki±  h;(^)  . 


(39)* 


Theorem  9.* 


j(N),  .  .  !0^ 

K  jur 


U) 

ilM 

(z)  ’ 

h  (z) 
n 

(N)*/  N 

m,  ,  -^-v  (-) 


/  m-1+6  „ .  m+6 

-  -  -n'^ - , -n+ - 

2  J  2  2 


=  (z/2) 


\-2n-l-v+6  ,m+v+l,-2n-l+6'*Tnl 

T  j-i:  ...  '^/A\ 


QPi”^-2n-l+6-v;-zV4) 


(N)*  (N) 

2  m,v(2)  = 


where  is  given  in  (45)*. 


Theorem  13.*  If  v  is  not  a  positive  integer  or  zero, 

.  2n+2  A 


•n+6  -v 


-z^/4 


m  .  2n+2-6+v 
=  sin  wr(2L^*v) 

r(2n+2-6+v)  '^2n+2-6+v^  ^  ^ 

(n/2)(z/2)^"''^'^^^ 

■*■  tan  vrr  r(2n+2-6+v)  '^m+v(2)J2ttf2-6+v(^)  • 

Equation  (60)*  can  te  rearranged  so  that  with  the  aid  of  L’Hospital’s 
theorem,  we  can  get  a  representation  of  the  error  when  v  becomes  a  posi¬ 
tive  integer  or  zero.  This  result  is  emitted.  Hovever,  for  arbitrary  v  , 
we  always  have 
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oFi(.2n-l+6-v;-274)E;  '  (z) 


,  ,  .2n+2  ,  ^ 

(z/2)  r(n+l-5+v) 

(n+l):r(2nt-2-6+v) 


(^)lF| 


1 

n+2,-n+6-v 


^(72) 


2n+2-6+v 


r(2n+2-6+N)) 


’^2n+2-6  +v  ( ^  ^  ) 


)((^/2)^1 


r(2n+2-6+v)r(2n+3-6+v) 


(61V 


■where  s  =  n-5+v  (s  =  »)  if  m  is  (is  not)  a  positive  integer  or  zero. 
Clearly  the  backward  recurrence  scheme  is  convergent.  Further,  for  n 
sufficiently  large,  n  »  m  ,  the  relative  error  is  essentially  indei>endent 
of  m  .  For  convenience  in  the  applications  we  record  the  formula 


m,v 


(z) 


(z/2)^^^7(n^l-6-Hv) 

(n+l):r(2n+2-6+v) 


r(2n+2-6+v)  7n-<-2-6+v^^^^in+v^^)  ' 


(62)# 


Remark ;  Let  v  ,  n  and 
m  only.  Then  £^^.^*(2) 
see  (21)#. 


(N)* 

z  be  fixed  so  that  Eni,v  (z)  is  a  function  of 
satisfies  the  recurrence  formula  for  > 


Let 


^^,1* =  Jm+v(z)  -  .i^v(z)  > 


J  (z>p(N)^(z;) 

'Po,v 


where  is  given  by  (37)*. 


(79)# 


(80)# 


3l» 


Bieorem  16.* 


(lO*.  .  2nt2-6+v'‘''^0,v 

°«.v  (“> - VF77 

fo,^  ("> 


s2n<-2-m-6 


_  r(m+v)(z/2)"“  .  /  X.*  /  s 

r(2n+2i+v)  '^2n+2-6+v^=')^,n(2)  > 


Im-i]  /_  m-1  m-2 

[  2  j[  2  2 

*  2^?  \l-v-m,v+l.l-i 

^  r  6-11/  J-1  ^6 

n-^  —  -n+ - ,  -n+- 

2  J  2  2 


F 
2  3 


-2n-l+6 -V ,v+l, -2n-l+6 


2\ 


,  0  <  m  s  2n<-l-6  .  (81)-i 


(  n)*  * 

Further,  ^  (z)  satisfies  the  recurrence  formula  for  tpjjj+^(z)  ,  see  (21)-> 
Thus 


oit.c.) = 4?i:v(^)  -  e*(^) . 


,(N)*. 

^0,v 


GA"/.  (z) 


0  ,  Oif^z) 


'^2n+2-6+v^’^^ 

(N)*,  . 

cp^  (z) 


(82)* 


Finally,  for  convenience  in  the  applications,  we  record  the  formulas 


/  /„v4n+4-26+2v-m  ,  . 

^(N)*,  ,  (z/2)  ^  r(m+v) 

^m,v  vz;  -  r(2n+2.6+v)r{2n+3-6+v)r(v+l) 


I¥l( 

/  m-1  m-2 

2  ’  2 

3  ' 

\^l-v-ra,v+l,l-m 

i 

Jv(z) 


HO(n"^)], 

(83)* 


(84)* 
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NIMERICAL  EXAMPLES 


Let 


N=5,  n=3,  6  =  1,2  =  2/3,  \>  =  1/3. 
Values  of  >  9^^^(z)  and 


'  ^ *7  \  OT««:k  tT'tir^n  ■?  r»  ■f.fi'hlo 


m 


6 

5 

4 

3 

2 

1 

0 


1 
19 
305 
3984 
40145 
2  84999 
11  80141 


9^^^(z)  =  680  75120/81 
n^^^(z)  =  1  38952  97360/6561 


Since 

•4^^=  0.77645  82114,  =  0.51341  71190  . 

r(v+i) 


the  Case  I  and  Case  II  approximations  are 


0.84272  08930  and  0.84272  10326  , 
respectively.  To  10  decimals 


11/2(2/3)  =  0.842772  08819  . 

Thus  the  errors  in  the  Case  I  and  Case  II  approximations  are 

-7  6 

-O.lll'lO  and  -0.15iao'  , 

respectively.  Using  (62)  aind  (76)  each  with  0(n“^)  and  the  term  involving 
K^^(z)  neglected,  the  approximate  Case  I  and  Case  II  errors  are 


and  -O.ug'lO"^  , 


respectively. 

For  a  second  example,  let 


N=5  ,n=3,6  =  l,  z  =  2  ,v=0 


Again  we  illustrate  the  Case  I  and  Case  II  schemes, 
data. 


m 


We  have  the  following 


5  6  =  611 

4  31 

3  130  =  4515 

2  421 


1 

972 

e2  =  7 

.38905  6099 

0 

1393 

i(N) 

m 

(?) 

z) 

m 

Case  I 

Case 

II 

0 

2.27986 

9067 

2.27972 

4285 

2.27958 

5302 

1 

1.59083 

4697 

1.59073 

3672 

1.59065 

6855 

n 

c. 

0.68903 

4370 

0.68899 

0613 

0.68694 

8448 

3 

0.21276 

5957 

0,21275 

2446 

0.21273 

9959 

4 

0.50736 

4975-10"^ 

0.50733 

2755-10“^ 

0.50728 

5700-10"^ 

5 

0.98199 

6727- 10"^ 

0.98193 

4365* 10'^ 

0.98256 

7932-10-2 

6 

0.163GC 

6121- 10”^ 

0.16365 

5728 '10-2 

0.16001 

7336-10-2 

Error 

Relative 

Error 

m 

Case  I 

Case  II 

Case  I 

Case  II 

0 

-0.284-10-^ 

-0.139-10~5 

-0.124-10-5 

-0.610-10-5 

1 

-0. 198- 10-^ 

-0.968-10"^ 

-0.124-10-5 

-0.609-10-^ 

-0.859-10-^ 

-0,422-10-^ 

-0.125-10-5 

-0.612- 10-*^ 

•z. 

-0.260.10-“^ 

-0.125  •10*'^ 

-0.122-10-5 

-0.587-10- " 

4 

-0.793-10-5 

-0.471-10-5 

-0.156-10-5 

-0.929-10-“^ 

5 

0.571-10“5 

0.634-10-5 

0.581-10-5 

0.645-10-^ 

6 

-0.365-10-^ 

-0.364-10-^ 

-0.228-10-1 

-0.228-10-1 
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> 

Here  the  entries  in  the  Ini(z)  column  are  correct  for  the  number  of  < 

decimals  given. 

Using  the  first  lines  of  (62)  and  (76),  each  with  0(n‘’^)  «md 
the  term  involving  lQn+v(z)  neglected,  the  approximate  relative  error 
for  Cases  I  and  II,  respectively,  are  -0.116*10“^  and  -0.537-10“^,  respec¬ 
tively. 

In  the  table  below,  we  record  the  approximate  errors  obtained  by 
use  of  (62)  with  0(n"^)  omitted  for  m  =  6  and  5  and  by  use  of  (21),  see 
the  remark  following  Theorem  13,  for  the  lower  values  of  m  .  This  is 
called  Case  I,  (62) -(21)  in  the  table.  We  also  present  the  analogous  Case  II, 

(76)-(21)  data.  In  each  instance  known  tabular  values  of  1^(2)  and  Iy(2) 
were  used.  In  practice,  we  suggest  using  (34)  or  the  lead  term  of  the 
uniform  asymptotic  expansion  of  {(^+.^(2)  developed  by  Olver  [9]  .  For 

l2nf2-6+\>(^)  ’  term  in  the  uniform  asymptotic  expansion 

for  this  function  which  is  also  given  in  the  source  just  cited.  We  also 
suggest  that  computation  of  the  gamma  functions  be  simplified  as  follows. 

With  |R(af)l  <  1  and  r  a  positive  integer,  we  have 


r(i^+l)  =  r: 


r(r4Ty^-l) 

r(r+i) 


rlr“[l+0(r-^)] 


and  for  r  sufficiently  large,  we  neglect  0(r"^)  .  The  approximation  is 
of  course  superfluous  if'y  =  0.  If  ct  =  ±  ^  ,  the  approximation  may  still 
be  used  though  known  tables  of  the  gamma  function  for  half  an  odd  integer 
may  be  preferred  lio].  If  more  precise  values  of  the  gamma  functions  are 
required,  see  [ll]  . 


Approximate 

Error 

m 

case  I,  (62)-(2l) 

case  II,  (76)-(21) 

0 

-0.264-10"^ 

-0.114.10-5 

1 

-0.184. 10’^ 

-0.765*10-^ 

2 

-0,797 -lO'^ 

-0.373.10-^ 

3 

-0.242. 10"^ 

-0.110-10-^ 

4 

-0.723.10"^ 

-0.410.10-5 

5 

0.475.10'^ 

0.536-10-5 

6 

-G. 310. 10"^ 

-0.309.10-^ 

For  a  final  example,  we  illustrate  Case  III  using  the  data  of  our 
second  example.  We  get  the  following  numbers. 
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m 

m 

m+v 

Approximate  Error 

(82)-(85)  _ W _ 

0 

2.27958 

5302 

0 

0 

0 

1 

1.59063 

6693 

0.162-10"® 

0.121-10"® 

0.276-10"® 

2 

0.68894 

8609 

-0.162-10"® 

-0.121-10"® 

-0.276-10"® 

3 

0.21273 

9475 

0,484-10"® 

0.363-10"® 

0.551-10"® 

4 

0.50730 

1826-10'^ 

-0.161-10"® 

-0.121-10"® 

-0.165-10"® 

5 

0.98187 

4502-10“^ 

0.693-10"® 

0,520-10"® 

0.661-10"® 

6 

0.16364 

5750-10-^ 

-0.363-10"^ 

-0.272-10"^ 

-0.331-10"'^ 

In 

the  above, 

the  approximate 

error  (82)-(83) 

means  that  G^+|(z)  is 

approximated  by  (83)  with  m  = 

=  1  and  0(n"“) 

neglected,  and  subsequent 

approximate  values  of  the  error  are  found  by  use  of  the  recursion  formula 
in  (82).  Use  of  the  recurrence  formula  in  this  fashion  is  stable  as  the 
magnitude  of  the  error  is  an  increeising  function  in  m  .  Also  Eq.  (84) 
means  this  equation  with  0(m“^)  and  0(n“^)  neglected. 

A  measure  of  the  accuracy  of  th:  three  schemes  treated  can  be  had 
by  use  of  normalization  relations.  Thus  if  the  Case  III  procedure  is  employed, 

(N) 

then  (38)  and  (46)  with  Ijj+^(z)  replaced  by  ijj+^(z)  are  available  as 

checks.  Similarly,  equations  (46)  and  (38)  are  available  as  checks  for  the 
Case  1  and  Case  II  techniques,  respectively.  For  some  other  useful  normali¬ 
zation  relations,  see  2,  pp.  45,  463. 

Analyses  of  the  error  in  the  backward  recursion  process  for  the 
solution  of  a  general  second  and  higher  order  linear  difference  equation 
have  been  given  by  a  number  of  authors.  Some  authors  have  studied  the 
ceise  of  Bessel  functions  directly.  We  make  no  attempt  to  survey  the  various 
contributions  here.  Pertinent  references  are  given  by  Wimp  [4I  .  Suffice  it 
to  say,  none  of  the  analyses  have  the  precision  and  simplicity  of  those 
developed  in  the  present  paper.  We  deliberately  chose  N  and  as  a  conse¬ 
quence  n  small  (W=5,n=3)  in  our  numerical  examples  to  put  our  asymptotic 
estimates  under  a  severe  test.  The  efficiency  and  realism  of  oxir  error 
formulas  is  manifest. 


CONCLUDING  REMARKS 


It  appears  that  the  techniques  developed  here  for  the  Bessel  func¬ 
tion  Iy(z)  can  be  extended  to  analyze  more  general  second  and  higher  order 
difference  equations.  In  particuleir,  it  would  be  useful  to  have  analogous 
results  for  confluent  forms.  This  we  intend  to  do 

in  future  papers. 


J 


I 

■ 
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